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Knowledge of the energy spectrum of *B neutrinos is an important ingredient for interpreting 
experiments that detect energetic neutrinos from the Sun. The neutrino spectrum deviates from the 
allowed approximation because of the broad alpha-unstable *Be final state and recoil order correc- 
tions to the beta decay. We have measured the total energy of the alpha particles emitted following 
the beta decay of B. The measured spectrum is inconsistent with some previous measurements, in 
particular with a recent experiment of comparable precision. The beta decay strength function for 
the transition from *B to the accessible excitation energies in *Be is fit to the alpha energy spec- 
trum using the R-matrix approach. Both the positron and neutrino energy spectra, corrected for 
recoil order efi^ects, are constructed from the strength function. The positron spectrum is in good 
agreement with a previous direct measurement. The neutrino spectrum disagrees with previous 
experiments, particularly for neutrino energies above 12 MeV. 

PACS numbers: 23.40.Bw, 23.60.+C, 26.65.+t, 27.20.-|-n 



I. INTRODUCTION 

The most carefully studied component of the solar 
neutrino flux is due to neutrinos from the /3^ decay 
of *B. The *B neutrinos account for most of the sig- 
nal in the Homestake ^^Cl neutrino capture experiment 
[l| and nearly all of the solar neutrino events in the 
Kamiokande |3], Super-Kamiokande Q, and the Sud- 
bury Neutrino Observatory (SNO) 3 water-Cherenkov 
experiments. Results from the SNO heavy water detector 
demonstrate the existence of a Vf^^r component of the so- 
lar neutrino flux jj| . The solar neutrino data is explained 
by flavor oscillations and non-zero neutrino mass 5] . The 
recent results of the KamL AND reactor Fg disappearance 
experiment |(| support the oscillation interpretation [5J- 
The neutrino oscillation solution implies that the solar 
®B Vf, energy spectrum is distorted. Knowledge of the 
primary ®B neutrino spectrum is a necessary ingredient 
for the proper interpretation of the solar neutrino data. 

A diagram illustrating the ®B decay chain is shown 
in Fig. n The ®B ground state (J'^=2"'") undergoes an 
allowed /J"*" transition to a broad range of excitation en- 
ergies in the a unstable ^Be daughter. We define the 
®B Z?"*" decay strength function as the probability that a 
given excitation energy in ^Be will be populated. The 
strength function is determined by measurements of the 
a particle energy spectrum following the breakup of the 
daughter ^Be nucleus, and is necessary to construct the 
neutrino spectrum. 

Transitions from *B to the ^Be ground state (J'^=0+) 
or the broad state at 11.4 MeV (J'^=4+) are second 
forbidden and highly suppressed. States in *Be with 
J'^=l+,3''" are not energetically accessible. The ®B /?+ 
decay thus proceeds exclusively through the resonance 
2+ structure in ®Be, described in the R-matrix formalism 




la, 

-0.092 MeV^g'O * MeV 

'He+'He 'Be 



FIG. 1: Nuclear levels in A=8 which lie below the *B ground 
state. Spins, parities, and excitation energies relative to the 



Be ground state are indicated. 



as a series of interfering 2+ states. We have measured the 
total energy of the a particles emitted following *B /?+ 
decay in a recent experiment 7], reviewed in Sec. ^] The 
data are analyzed in the framework of the many-level R- 
matrix approximation to fit the /3+ decay strength func- 
tion, which is presented in a table and compared with the 
results of the previous precision measurement [31 in Sec. 
mil The R-matrix approach employed here is not essen- 
tial for obtaining the neutrino spectrum, but provides a 
convenient way to characterize the experimental data. 

The *B neutrino spectrum is subject to corrections due 
to recoil order matrix elements. Measurements involving 



the f3~^{f3~) decay of ^B(^Li) and the radiative decays 
of the 2+ doublet in ^Bc with excitation energies near 
16 MeV, shown in Fig. Q] are used to extract the recoil 
order matrix elements which contribute to ®B /?+ decay. 
A review of past recoil order measurements is presented 
in Sec. |E| 

Both the ®B positron and neutrino energy spectra are 
deduced, using the strength function and applying recoil 
order and radiative corrections, in Sec. The agree- 
ment between this work and the previous direct mea- 
surement of the positron spectrum 9] is discussed, and 
the neutrino spectrum is presented in a table. 
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II. 



THE ALPHA SPECTRUM MEASUREMENT 



A description of the a-spectrum measurement dis- 
cussed here has been presented previously |3- In this 
section we briefly recount the experimental technique, 
focusing on the experimental uncertainties. 



A. Experimental Technique 

A beam of ^B ions was implanted near the midplane of 
a 91 /im thick planar Si detector. An implanted source 
eliminates the possibility of energy loss outside the sensi- 
tive region of the detector, a systematic effect in all pre- 
vious experiments. The detector thickness was just suffi- 
cient to stop a particles emitted with the highest possible 
energy (about 8.5 MeV). Thus the full energy of both a 
particles was detected while the positrons, usually close 
to minimum ionizing, deposited only a small amount of 
energy. The systematic effect of positron energy was fur- 
ther reduced with a coincidence detector, selecting events 
where the positron trajectories were close to normal to 
the Si detector surface. The system was calibrated using 
a beam of ^"Na which was also implanted near the detec- 
tor midplane. The ^°Na decays with 20% probability to 
a unstable levels in ^'^Ne, providing calibration lines of 
well-known energy jl3|. An external ^^^Th a source was 
used to provide additional calibration lines. 

The experiment used the ATLAS superconducting lin- 
ear accelerator at the Argonne National Laboratory. The 
®B (ti/2='''70±3 ms) beam was produced by the In-Flight 
Technique ill(| using the ■^He(^Li,*B)n reaction. The pri- 
mary ^Li beam, with energy 36.4 MeV, bombarded a 3.5 
cm long gas cell filled with 700 mbar ^He and cooled to 82 
K. The gas cell was separated from the beam-line vacuum 
by titanium windows. The pressure and temperature in 
the cell were held constant to ±1%. 

Fully stripped *B products were separated from the 
primary beam with a 22° bending magnet, and trans- 
ported through an Enge Split Pole spectrograph. A gas- 
filled detector located in the focal plane of the spectro- 
graph [l^ identified the ^B products by mass, nuclear 
charge, and energy. The spectrograph was then adjusted 




FIG. 2: The experimental setup used to produce the *B 
(^"Na) beam, separate it from the primary ^Li (^^F) beam, 
and select ions with energy 27.3 MeV (170 MeV) for implan- 
tation into the Si detector (not to scale). 



so that ^B ions with energies of 27.3±0.2 MeV were inci- 
dent on the planar Si detector with a 150 mm^ active area 
(13.8 mm diameter), located adjacent to the focal plane 
detector. An 11 mm diameter Ta collimator masked the 
edges of the detector. The beam was cycled (1.5 sec 
on/ 1.5 sec off) and data taken only during the beam-off 
cycles. The average implantation rate was 3 ^B ions/sec, 
and 4.5 X 10^ decays were observed over six days. 

Using the gas-filled focal plane detector, the *B beam 
purity was measured to be about 10""^. A portion of the 
low-energy tail of the primary ®Li beam, as well as a par- 
ticles, deuterons, and protons with the proper magnetic 
rigidity to traverse the spectrograph, accounted for the 
remainder of ions incident on the detector. A ^Li beam 
incident on ''He cannot produce any /3 delayed particle 
emitters other than ^B, which could create a background 
during the beam-off data acquisition cycles. No products 
resulting from possible interactions between ^Li and the 
titanium windows of the gas-cell were observed with the 
proper rigidity to be incident on the Si detector. 

The /3 particle detector, located 12 mm behind the 
Si detector, was a 25 mm diameter x 2 mm thick plastic 
scintillator coupled by a lightguide to a Hamamatsu R647 
photomultiplier tube. The detector identified a subset of 
events where the positron from the *B decay exited the Si 
detector with a trajectory within 30° to normal. Roughly 
16% of the observed events occurred in coincidence with 
a count in the (3 detector, consistent with expectations 
from detector geometry. The Si/scintillator detector sys- 
tem was cooled to -5°C. A schematic representation of 
the apparatus is shown in Fig. |21 

The calibration using implanted ^°Na was performed 
immediately before the ^B run. The ^"Na (3~^ delayed 
alpha particles provided three calibration lines near the 



region of the ^B a spectrum peak, with energy releases of 
2691.9±1.2, 3099.0±2.2, and 5544.0±2.8 keV O- The 
^°Na (ti/2=448±3 ms) beam was produced by using the 
i9p('3jjg^2n)^''Na reaction and separating fully stripped 
^•^Na ions with energies of 170.0±1.5 MeV. A mylar de- 
grader foil of thickness 85±4 /im in front of the detector 
slowed the ions prior to implantation. As in the ^B runs, 
the beam was cycled (1.0 sec on/ 1.0 sec off). An average 
implantation rate of 7 ^"Na ions/min was achieved, and 
over one day 1.0x10'* decays were observed. The raw 
energy spectra from the *B and ^"^Na runs are displayed 
in Fig. 

The integrated incident flux on the Si detector, mon- 
itored by the spectrograph focal plane detector, was an 
order of magnitude below threshold for detector damage 
|13J |. and no gain variation from damage was expected. 
The gain was monitored with the centroid of the ®B a 
spectrum and was found to fluctuate within ±0.25%, cor- 
responding to ±7 keV at the spectrum peak. The fluc- 
tuations are about two times larger than the statistics. 
External a particle sources were not reliable for moni- 
toring gain shifts because of the accretion of residual gas 
onto the cooled Si detector, degrading the a particle en- 
ergies by 10-20 keV over the course of the seven day run. 
The accreted gas was not sufficient to appreciably de- 
grade the ^B and ^°Na beams incident on the detector, 
and did not affect implantation depth. 

The electrostatic sweeper used to stop the beam was 
not perfectly efficient, allowing a weak beam during the 
counting cycles. Protons with energies near 8.7 MeV, 
produced in reactions from the primary beams, had the 
right rigidity to traverse the spectrograph and hit the 
Si detector. The 8.7 MeV protons passed through the 
Si detector and into the /? detector, producing a peak 
near 800 keV in the coincidence data. These protons 
were rejected based on the large pulses observed in the 
/3 detector, which were much larger than the pulses from 
minimum ionizing positrons. Any ions heavier than pro- 
tons with the proper rigidity to reach the Si detector were 
stopped |l4| and rejected by the coincidence requirement. 

The analysis provided here is slightly improved over 
Ref. Q. Proton events were removed by cuts, as noted 
above. Random coincidences were removed with cuts on 
the time spectrum recorded between the /3 detector start 
and Si detector stop signal. Additional cuts on the time 
spectrum eliminated a small amount of background from 
external /3 decay activity. 



B. Experimental Uncertainties 

The largest sources of experimental uncertainties in 
the a. spectrum measurement, affecting the ^B neutrino 
spectrum, are: (1) The temporal gain variation observed 
over the seven days of data collection. (2) The uncer- 
tainty in correcting the energy deposited by positrons, 
which includes the uncertainty in implantation depth of 
the ^B and ^"Na ions. (3) The uncertainty in the energy 
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FIG. 3: The measured *B /?^ delayed two alpha spectrum 
shown with the ^°Na /S"*" delayed alpha lines used for calibra- 
tion. The data shown here correspond to events coincident 
with the beta detector. 



scale calibration. 

As noted, the nonstatistical gain variation over the 
seven day run was of magnitude ±0.25% which corre- 
sponds to a ±7 keV uncertainty at the peak of the a 
spectrum. A gain correction was not applied, instead an 
uncertainty is included in the energy scale. This is the 
dominant source of uncertainty in the measurement. 

The depth and distribution of the implanted ions was 
estimated using the TRIM Monte Carlo simulation [Ij|. 
The uncertainty is taken as ±4.6%, the average devia- 
tion of TRIM estimates from measured stopping power 
[Ti| . For ^B ions of incident energy 27.3±0.2 MeV, TRIM 
predicts an average implantation depth is 42.2±2.0 /im. 
The full-width half maximum of the implantation depth 
is 0.7 ^m. For ^ONa ions of energy 170.0±1.5 MeV, 
first passing through the mylar degrading foil of thick- 
ness 85±4 fim, TRIM predicts an average implantation 
depth of 48±6 /im for ^°Na ions, with a full width half 
maximum of 1.3 /im. 

Uncertainties in ion implantation depth correspond to 
uncertainties in energy deposited by positrons. On av- 
erage, minimum ionizing positrons deposit 0.6 keV//im 
in Si, so that in the case of ®B (^°Na) the uncertainty in 
implantation depth corresponds to an energy uncertainty 
of ±1.2 keV (±3.6 keV). 

The positron energy loss in the Si detector was esti- 
mated using the EGSnrc simulation 15j. Simulations 
account for detector geometry and positron energy spec- 
tra and assume the ranges of ion implantation depths 
discussed previously. Probability distributions for energy 
loss by positrons were obtained for the subset of data as- 
sociated with a coincidence count in the (3 detector, and 
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FIG. 4: Fits to the ^''Na calibration lines. The open circles 
indicate the coincidence data set, while the solid circles indi- 
cate the total data set, i.e. no coincidence requirement. The 
curves show the best fit function, described in the text. 



for the total data set. The uncertainty associated with 
these simulations was estimated by comparing the total 
®B data set to the coincidence data set. The effect of the 
positron correction lowered the a-spectrum peak of the 
total data set by 55 keV, and the peak of the coincidence 
data set by 24 keV. After the correction, the peaks of the 
two data sets agreed to within 2 keV. The uncertainty as- 
sociated with the simulation is thus assigned as ±2 keV. 
The use of the total data set to estimate uncertainty in 
positron energy loss was not compromised by the beam 
leakage background, since beam particles with the proper 
rigidity to hit the detector had energies far from the ^B 
spectrum peak. 

The average pulse height defect of the recoil ^^O nuclei, 
which carry one fifth of the energy of the a disintegra- 
tions following ^'^Na decay, has been directly measured 
for -'^^O nuclei in the energy range of interest i()\. The 
correction is 40-50 keV for the various ^"Na alpha lines, 
with an uncertainty of ±5 keV. The TRIM Monte Carlo 
simulation ^^ was used to model the ionization energy 
loss of ^^O in silicon, and agreed within 2 keV with the 
average values of ionization loss observed in [l6J. We 
have applied the TRIM results, scaled by 2 keV to agree 
with the experimental results, to approximate the pulse 
height spectrum of ^^O nuclei in a silicon detector. 

The energy spectrum from the ^''Na decay was used to 
calibrate the energy scale. The calibration lines were fit 
to the pulse height spectrum predicted by TRIM, con- 
voluted with the positron energy loss distributions and a 
Gaussian component to approximate detector noise. The 
position and amplitude of the lines were free parameters, 
as well as the Gaussian width. Results of the fit to two 
of the lines, resulting from ^°Na /3+ decays which led to 
a energy releases of 2691.9±1.2 and 3099.0±2.2 keV, are 
shown in Fig. 21 for both the total and coincidence data 
sets. 

Pulser tests performed before and during the data col- 




5 6 7 

Energy (MeV) 



10 



FIG. 5: Residuals from the calibration process. The solid cir- 
cles show the residuals for the three ^'^Na calibration lines to 
the best linear fit (energy vs. ADC channel). The open circles 
correspond to the six external ^^^Th source alpha lines, which 
were not used in the calibration fit. The thin solid curves are 
the la error bands associated with the '^'^Na calibrations. The 
dashed curves are the la error bands of a separate calibra- 
tion from the ^^^Th source. The thick solid curves show the 
total la uncertainty in the energy scale, which is significantly 
larger than the calibration uncertainty alone due to temporal 
gain variation. 



lection period indicate a negligible quadratic component 
in the relationship between pulse height and ADC bin, so 
the ^"Na calibration was performed using a two param- 
eter linear fit to the three dominant lines along with the 
zero energy ADC channel precisely determined by pulser 
tests. The external ^^^Th source emitted a particles at 
six distinct energies, 5.341, 5.423, 5.686, 6.288, 6.779, 
and 8.784 MeV 10] and was used to perform an indepen- 
dent calibration. Data used for the ^^*Th calibration was 
taken immediately after the detector was placed in vac- 
uum and cooled, before an appreciable amount of resid- 
ual gas condensed on the detector surface. The a particle 
energies were corrected for energy loss in the source and 
the measured 27±4 /ig/cm^ detector dead layer. The 
magnitude of the corrections for the various lines was 31- 
38 keV, with a characteristic uncertainty of 4-5 keV. A 
comparison of the residuals from the two calibrations is 
shown in Fig. [S] The figure also shows the total uncer- 
tainty in the energy scale, dominated by the temporal 
gain variation. 



III. R-MATRIX REPRESENTATION OF THE 
ALPHA SPECTRUM 

The many-level R-matrix formalism has previously 
been used to parametrize data from nuclear processes 
involving ^Be, in particular the alpha spectrum following 
«B /?+ decay @,[l3,[ll[ll[23. The application of the R- 
matrix to /3 decay is an approximation, and the physical 



significance of R-matrix fit parameters is not clear. 

In principle, it is possible to deduce the ^B /9+ decay 
strength function directly from the measured a spectrum 
without resorting to R-matrix formalism. The R-matrix 
approach, however, gives a good fit to the observed a 
spectrum and provides a convenient method for propa- 
gating systematic uncertainties in the a spectrum to the 
neutrino spectrum. Systematic uncertainties in the a 
spectrum dominate the statistical uncertainties, justify- 
ing the representation of the data by a smooth function. 



A. Energy Levels in *Be 

In the R-matrix approximation, the j3~^ decay of ®B 
proceeds with varying strength through a region of inter- 
fering nuclear states in *Be which immediately decay into 
a particles. Each state is characterized by an excitation 
energy, £j, a reduced width, 7^, and a P~^ decay strength 
quantified by Fermi and Gamow- Teller matrix elements, 
Mpj and Mgtj- 

As discussed in Sec. Q] only 2+ states in ^Be are consid- 
ered. A numerically accurate R-matrix fit to the observed 
a spectrum requires the three 2"*" states in ®Be shown in 
Fig. n as well as one phenomenological background state 
approximating the combined effect of all higher-lying 2"*" 
states. It has been shown |1^^J that R-matrix fits using 
only these four states were sufficient to describe previous 
experimental data. We have repeated the analysis dis- 
cussed in Ref. ^3, which explicitly included a greater 
number of 2+ states, and have verified that the four state 
R-matrix approximation is sufficient to describe the a 
spectrum reported here. 

The state labeled (j=l), with excitation energy near 
3 MeV and width of about 1.5 MeV, is responsible for 
the peak of the observed a spectrum. The excitation 
energy, £1, and reduced width, 71, are considered free 
fit parameters. Shell- model considerations [2l| indicate 
no significant Fermi decay strength to this level, as dis- 
cussed in Ref. |22| which reports measurements of the 
(3-v-a correlations in ®B and ^Li consistent with a pure 
Gamow- Teller decay. We take the Fermi decay strength 
to vanish, Mfi=^, while the Gamow- Teller matrix ele- 
ment, M-GTi, is a free parameter. 

The next two states (j=2,3) form a nearly degen- 
erate doublet with excitation energies 16.626(3) and 
16.922(3) MeV [l3| which are well known to be almost 
maximally mixed in isospin. We describe the isospin mix- 
ing of the doublet using the standard formulation ^_17j and 
consider the energy eigenstates "02 and -03 in terms of the 
isospin eigenstates (J)a (T=0) and 0_b (T=1), 



i/'2 = a4>A + , 



03 = P4>A - a(t)B, 



(1) 



where a and /? are mixing parameters with a^ + f3'^ — 
1. Since a decays from a T=l state are forbidden, the 
parameters a and /? may be approximated from the level 
widths, 



An accurate R-matrix description of the alpha spec- 
trum requires a,/9 >0. The energies, £2 and £3, and 
reduced widths, 72 and 73, of the doublet are well con- 
strained by a-a scattering experiments |23l | and are held 
constant. 

The decomposition of the doublet into its component 
isospin eigenstates allows a simplified description of the 
Fermi and Gamow- Teller strengths. The T=0 state, (J>a, 
has a Gamow- Teller strength treated as a free parame- 
ter, Mgta- The T=l state, 0b, is the isospin analog of 
the ^B and ®Li ground states and is populated by Fermi 
decay with a strength given by the superallowed Fermi 
matrix element, M.fb=V^- The Gamow- Teller decay to 
the T=l component, or Fermi decay to the T=0 com- 
ponent, may be nonzero due to isospin breaking but has 
been estimated to be negligible ^J in this context, as 
discussed in Ref. [ig- Hence we take AicTB—O and 
■Mfa—0. The matrix elements of the isospin eigenstates 
are then related to the matrix elements of the energy 
eigenstates by Eq. ^ 
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(3) 



(4) 



The background state labeled (j=4) has an excitation 
energy held fixed to a value near that used in recent 
works |1, liallOl' f 4=37.0. The parameter £4 could be 
allowed to float, but the quality of the fit is very weakly 
dependent on its value. The reduced width, 74, and 
the Gamow- Teller matrix element, Mgta, are free pa- 
rameters. The Fermi strength is taken to be negligible, 
Mf4^0. 



B. Form of the R-matrix Function 

The R-matrix approach gives a parametrization of the 
*B /S"*" decay strength function, indicating the probability 
that ®Be is populated at a give n excitation energy, Ea,. 
The function takes the form Il3l 



dN 
dE^ 



Nt 



1/2 



6166 



UiE^)(^a^{E^)+c\E^)y (5) 



Here N is the total number of observed decays. 



N = 



dN 
dE. 



-dE^ 



(6) 



a 



r2/(r2 + r3), /^^ = r3/(r2 + r3). 



(2) 



The lifetime of ^B, ti/2, is 770±3 msec l^- The unit- 
less integrated phase space available to the (3 decay lep- 
tons, /^(Ea;), including the Fermi function and outer ra- 
diative corrections, has been evaluated according to the 
parametrization given by Wilkinson and Macefield |25| . 
The Fermi and Gamow- Teller matrix elements, a(Ea;) and 
c(E2;), are parametrized by 
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FIG. 6: Top panel: R-matrix fit to the observed decay spec- 
trum. Bottom panel: Residuals to the fit, scaled by the square 
root of the fit value. 



The P(E2;) and S(Ea;) are the penetrability and shift 
factor arising from the regular and irregular solutions of 
the Coulomb equation of L=2 a particles, defined in Ref. 
|2g . As in previous works H, HI, 123 , we evaluate the 
Coulomb functions at matching radius rc=4.5 fm, and 
choose the boundary condition, B, to be S(£i). 



C. Application to Data Set 

In cases where the *B decays at rest, the recoil of the 
daughter ^Be nucleus will cause the total energy spec- 
trum of the emitted a particles to deviate from the f3^ 
decay strength function given in Eq. For a given ex- 
citation energy of ^Be, the recoil energy distribution is 
exactly calculable and takes an average value of 7 keV at 
the most probable excitation energy near 3.0 MeV. 

In addition to accounting for the ^Be recoil, the 




Be Excitation Energy (MeV) 



FIG. 7: The B /3 decay strength function, determined by 
R-matrix fits to the *B alpha spectrum presented in this work 
(black) and in Ortiz et al. |a] (grey). The spectra are scaled to 
share the same peak height. The insert shows the locations 
of the spectrum peaks, on which the neutrino spectrum is 
highly dependent. The width of the lines in the insert indicate 
the magnitude of the ±1(t experimental uncertainties. The 
thin feature in the blue curve arises because the dominant 
uncertainty is a multiplicative factor in the energy scale. 



strength function (Eq. [SJ must be convoluted with 
the probability distribution of energies deposited by the 
positron, discussed in Sec. ^ The detector line shape, 
approximated as a Gaussian with width 25 keV, deter- 
mined by fits to the ^°Na data sets, was also included but 
had a negligible impact on the fit due to the large width 
of the a spectrum. 

The a spectrum data was fit using the log- 
likelihood minimization function [j^l- The best fit gave 
X^/dof=3249. 7/3376, indicating a satisfactory fit. The 
best fit parameters are given in Table HJ and the best fit 
is compared to the data in Fig. El The strength function 
is presented in numerical form in Table Ull 

The strength function reported here disagrees with the 
result of Ortiz et al. ,8|. For both measurements, the un- 



TABLE I: Values of R-matrix parameters determined by a fit 
to tlie coincidence a spectrum data, using a matching radius 
of rc=4.5 fm. 



Parameter 



Value 



£i 

£2 

£3 

£4 

7? 

7I 

7I 

7I 

Man 

Mgta 

Mgta 

Mfb 

Mgtb, Mfi, Mfa, Mf4 



3.043 MeV 

16.626 MeV" 

16.922 MeV" 

37.0 MeV 

1.087 MeV 

10.96 keV" 

7.42 keV 

5.619 MeV 

-0.1462 

2.423 

-0.1320 

0" 



"These parameters were held constant during the fitting proce- 
dure. 



certainty in the inferred neutrino spectrum is dominated 
by systematic effects. Smooth R-matrix fits to the alpha 
spectra thus provide a convenient way to compare the 
two results. Fig. [7| shows a comparison of the present 
results and a fit to the data of Ortiz et al. Q- Uncer- 
tainties in the Ortiz et al. curve are taken directly from 
Ref. li. 



D. Propagation of Systematic Uncertainties 

The R-matrix approach was used to propagate the sys- 
tematic uncertainties in the a spectrum measurement to 
the neutrino spectrum. As discussed in Sec. ^ the la 
uncertainty in the energy scale is characterized by a mul- 
tiplicative factor of 1± (0.275%), corresponding to about 
8 keV at the spectrum peak, added in quadrature with 
a constant offset of 3 keV. R-matrix fits were performed 
to the a spectrum using the iblcr energy scales, and the 
resulting zblcr strength functions were used to produce 
±lf7 neutrino spectra. 

An additional uncertainty was imposed to account for 
the rapid drop off of the a spectrum at low energies, 
where statistics are not sufficient to determine the spec- 
trum shape. The penetrability factor, P(Eix), is respon- 
sible for the drop off. The best R-matrix fit used P(Ej;) 
calculated for a matching radius of 4.5 fm. We approxi- 
mate the uncertainty at low energies by calculating P(Ej;) 
at 4.0 and 5.0 fm, the ±1ct matching radii recommended 
in Ref. [13 j and perform fits under these conditions. 
There is a strong dependence [13 between the energy 
of the background state, £4 and matching radius, re, so 
the parameter £4 was allowed to fioat for these fits. We 
note that the variation of matching radius is a signifi- 



cant source of uncertainty only for neutrinos at very high 
(E,, >15 MeV) and low (Ei. <0.5 MeV) energies. 



IV. RECOIL ORDER CORRECTIONS TO THE 
NEUTRINO SPECTRUM 

A. Background 

A proper description of *B f3~^ decay includes recoil 
order effects which cause, for example, the energy spec- 
tra and angular correlations of decay particles to deviate 
from the allowed approximation. Deviations are of order 
Eo/m„, where Eq is the positron cndpoint energy and 
m„ is the nucleon mass. The '^B /3+ decay has a particu- 
larly large endpoint energy (most probable Eq «13 MeV) 
and a small Gamow- Teller strength (log ft = 5.6) for an 
allowed decay. Recoil order effects in ^B are thus large 
compared to other nuclear systems. 

Measurements of the radiative decay of the ^B isospin 
analog state in ^Be [H IH S EH, and of the angu- 
lar correlation between and a particles emitted in the 
decays of ®B and *Li |23, HI, HI], determine the recoil 
order matrix elements. These results were first explicitly 
applied to the neutrino spectrum in 9, 34], where they 
were found to contribute at the 5% level. A more re- 
cent determination of the neutrino spectrum by Bahcall 
et al. i35il emploved the same recoil order treatment as 
in Ref. [33|. Bahcall et al. provided a conservative esti- 
mate of the uncertainty associated with the recoil order 
correction on the neutrino spectrum by setting the 3a 
uncertainty equal to size of the correction itself. A more 
recent determination of the neutrino spectrum by Ortiz 
et al. Q, applied recoil order corrections very similar to 
those in Ref. |3i |. 

The two most recent a spectrum measurements, by 
Ortiz et al. 8] and the one reported here, involved de- 
terminations of the energy scale significantly more precise 
than the measurements used by Bahcall et al. 35] . Also, 
a recent precision measurement of radiative decay in ®Be 
|3l| provides additional information on recoil order ef- 
fects, but has not yet been applied to the *B neutrino 
spectrum. In light of these recent experiments, recoil or- 
der effects are considered here with careful attention to 
the assignment of realistic uncertainties. 

The Fermi matrix element plays a small role in the 
(3'^ decay of ^B, contributing only to decays proceed- 
ing through the highest excitation energies in ^Be, as 
explicitly discussed in Ref. [l^. These low energy /3''" 
decays have suppressed recoil order corrections, and pro- 
duce neutrinos which have no impact on solar neutrino 
experiments. Consideration of the Fermi matrix element 
is thus omitted. 

A model independent treatment of recoil order ef- 
fects is given by Holstein |3y], whose notation we adopt 
here. Matrix elements contributing to the /3 decays of 
^B and ^Li are denoted by c (Gamow- Teller), b (weak 
magnetism), d (induced tensor), /, g (vector second- 



TABLE II: The *B /3^ decay strength function, as determined by fitting the experimental a spectrum to Eq. |3 The strength 
function is normahzed to 1000 when integrated with respect to MeV. Note that the energy spacing of data points varies to 
allow a more detailed description of the strength function near the peak. Uncertainties in the function are well approximated 
by deforming the energy scale. The \a uncertainty in the energy scale is characterized by a multiplicative factor of 1± (0.275%), 
corresponding to about 8 keV at the spectrum peak, added in quadrature with a constant offset of 3 keV. 



E. 


dN/dE^ 


E, 


dN/dE^ 


E. 


dN/dE^ 


E. 


dN/dE^ 


E, 


dN/dE^ 


E, 


dN/dE^ 


E, 


dN/dE^ 


E. 


dN/dE^ 


0.00 


0.000 


1.90 


67.991 


2.60 


362.922 


3.06 


462.243 


3.52 


315.256 


4.70 


117.714 


8.00 


24.113 


12.60 


2.349 


0.10 


0.000 


1.95 


78.048 


2.62 


374.309 


3.08 


457.936 


3.54 


309.037 


4.80 


110.362 


8.20 


22.147 


12.80 


2.054 


0.20 


0.004 


2.00 


89.464 


2.64 


385.390 


3.10 


453.156 


3.56 


302.941 


4.90 


103.693 


8.40 


20.333 


13.00 


1.788 


0.30 


0.024 


2.05 


102.397 


2.66 


396.095 


3.12 


447.953 


3.58 


296.971 


5.00 


97.622 


8.60 


18.655 


13.20 


1.548 


0.40 


0.081 


2.10 


117.004 


2.68 


406.353 


3.14 


442.380 


3.60 


291.129 


5.10 


92.074 


8.80 


17.103 


13.40 


1.333 


0.50 


0.207 


2.15 


133.440 


2.70 


416.096 


3.16 


436.485 


3.62 


285.417 


5.20 


86.986 


9.00 


15.667 


13.60 


1.141 


0.60 


0.440 


2.20 


151.844 


2.72 


425.257 


3.18 


430.314 


3.64 


279.834 


5.30 


82.305 


9.20 


14.334 


13.80 


0.970 


0.70 


0.832 


2.25 


172.320 


2.74 


433.776 


3.20 


423.911 


3.66 


274.380 


5.40 


77.984 


9.40 


13.100 


14.00 


0.818 


0.80 


1.443 


2.30 


194.916 


2.76 


441.597 


3.22 


417.319 


3.68 


269.056 


5.50 


73.984 


9.60 


11.956 


14.20 


0.685 


0.90 


2.348 


2.32 


204.543 


2.78 


448.671 


3.24 


410.576 


3.70 


263.859 


5.60 


70.270 


9.80 


10.896 


14.40 


0.568 


1.00 


3.640 


2.34 


214.498 


2.80 


454.958 


3.26 


403.717 


3.72 


258.789 


5.70 


66.814 


10.00 


9.915 


14.60 


0.466 


1.10 


5.435 


2.36 


224.770 


2.82 


460.423 


3.28 


396.777 


3.74 


253.845 


5.80 


63.588 


10.20 


9.006 


14.80 


0.378 


1.20 


7.878 


2.38 


235.344 


2.84 


465.045 


3.30 


389.783 


3.76 


249.024 


5.90 


60.572 


10.40 


8.166 


15.00 


0.302 


1.30 


11.155 


2.40 


246.203 


2.86 


468.808 


3.32 


382.763 


3.78 


244.323 


6.00 


57.744 


10.60 


7.389 


15.20 


0.238 


1.40 


15.507 


2.42 


257.323 


2.88 


471.708 


3.34 


375.742 


3.80 


239.742 


6.20 


52.589 


10.80 


6.672 


15.40 


0.185 


1.50 


21.242 


2.44 


268.676 


2.90 


473.750 


3.36 


368.740 


3.90 


218.531 


6.40 


48.008 


11.00 


6.011 


15.60 


0.140 


1.55 


24.750 


2.46 


280.230 


2.92 


474.946 


3.38 


361.778 


4.00 


199.905 


6.60 


43.909 


11.20 


5.402 


15.80 


0.104 


1.60 


28.761 


2.48 


291.947 


2.94 


475.318 


3.40 


354.870 


4.10 


183.532 


6.80 


40.221 


11.40 


4.842 


16.00 


0.076 


1.65 


33.347 


2.50 


303.785 


2.96 


474.895 


3.42 


348.031 


4.20 


169.102 


7.00 


36.886 


11.60 


4.328 


16.20 


0.056 


1.70 


38.584 


2.52 


315.694 


2.98 


473.712 


3.44 


341.275 


4.30 


156.340 


7.20 


33.856 


11.80 


3.857 


16.40 


0.044 


1.75 


44.561 


2.54 


327.623 


3.00 


471.809 


3.46 


334.611 


4.40 


145.009 


7.40 


31.093 


12.00 


3.426 


16.60 


0.097 


1.80 


51.378 


2.56 


339.512 


3.02 


469.230 


3.48 


328.049 


4.50 


134.904 


7.60 


28.566 


12.20 


3.033 


16.80 


0.001 


1.85 


59.146 


2.58 


351.300 


3.04 


466.025 


3.50 


321.595 


4.60 


125.854 


7.80 


26.246 


12.40 


2.674 


17.00 


0.000 



forbidden), j2, j3 (axial second- forbidden) , and h (in- 
duced pseudoscalar) . Since the decays proceed to the 
broad continuum in ^Be, the matrix elements should 
be considered as functions of the ®Be excitation energy, 
E.T,. Previous determinations of the ^B neutrino spectrum 
|3,l2Jjl23 neglected this energy dependence. 



B. Beta and Neutrino Energy Spectra 

The positron energy spectrum from an allowed de- 
cay proceeding between two energetically sharp nuclear 
states is given by 






■ ppEpiEo - EfifF{~Z, Ef3)R{Ef3, Eo)C{Ef3, Eo). 

(9) 



Here p^ and E^ are the momentum and total energy 
of the positron, and Eq is the positron endpoint energy. 
F(-Z,E^) is the Fermi function, which depends on the 
charge, Z, of the daughter nucleus and is negative for 
positron decays. The radiative corrections are contained 
in R(E^,Eo), which will be discussed in Sec. The 
recoil order effects are contained in C(E^,Eo), which has 
the form 



CiEp,Eo)^l 



2Eo 
SArrin 



d 
c 



2Ep 
SArrin 



3AmnEa 



d 
c 



hEo 
c 



En 



2Amn 



(10) 



where A=8 is the mass number. In the case of ^B the 
recoil order matrix elements are dependent on the ®Be 
excitation energy E^; (Ea;=A-Eo), where A=17.468 IVIeV 
is the total energy released in the ®B /3-a decay chain. 
(This discussion of positron and neutrino energy spectra 
ignores, for the sake of simplicity, the kinetic recoil of the 
daughter nucleus. This effect is included in the numer- 
ical calculations.) The positron spectrum is calculated 
by integrating Eq. |2]over all excitation energies in ^Be, 
weighted by the strength function determined in Sec. 11111 
The neutrino spectrum is obtained by the simple substi- 
tution Ei^=Eo-E^, and the application of different radia- 
tive corrections, discussed in Sec. 



C. Radiative Decay Measurements in Be 

The weak magnetism matrix element, 6, exerts the 
greatest influence on the neutrino energy spectrum. Its 
value is best determined under the strong conserved vec- 
tor current (CVC) hypothesis by measurements of the 
radiative decays of the *B isospin analog state in ^Be 
which, as discussed in Sec. IIIII is mixed between the two 
states of an energy doublet. The radiative decay is shown 
schematically in Fig. ^ 

In previous experiments [23, 123, Hfl 113 , a ^He beam 
was directed at a ''He gas cell to excite the doublet 
in '^Be. The ''He(*He,7)^Be cross section was mea- 
sured as a function of beam energy and angle of emis- 
sion of the 7 ray. These measurements determine the 
widths of the isovector Ml and E2 transitions, Fj^^ and 
i5i=F^^-'/F^-', as well as the widths of the isoscalar 
transitions, e=Tlf^° /Tlf^^ and 5o^TI=° /rlf^\ The ra- 
diative widths are considered as functions of Ej, . 

CVC relates the isovector radiative widths in *Be to 
the vector recoil matrix elements contributing to *B /3 
decay, 6, /, and g, 



A summary of the experimental results is given in Ta- 
ble mil The experimental results for the isoscalar con- 
tributions to the decay, e and Sq, agree with each other 
and are of the same order as various shell model pre- 
dictions compiled in Ref. |3l| . The experimental values 
for (5i from two of the experiments l2i^ |31l are in agree- 
ment, but differ from the results in |8nl | by about 3cr. 
The present work will use the more recent value of 6i 

TABLE III: Experimental determinations of the isovector and 
isoscalar Ml and E2 transition strengths. All quantities listed 
are integrated over final state excitation energies in *Be. 



Observable 


Experiment 


Value 


5i 


Ref. [28] (1975) 
Ref. [30] (1978) 
Ref. [31] (1995) 


0.045±0.Q27 

0.14±0.03'' 

0.01±0.03 


So 


Ref. [30] (1978) 
Ref. [31] (1995) 


0.26±0.03" 
0.22±0.04 


e 


Ref. [30] (1978) 
Ref. [31] (1995) 


0.00±0.03" 
0.04±0.02 


-l- Ml 


Ref. [29] (1977) 
Ref. [30] (1978) 
Ref. [^ (1996) 


4.1±0.6 eV' 

3.6±0.3 eV' 

2.80±0.18 eV 



b{E,) = Amr^^6rl=\E,)/{aE^^), (11) 



"These values are based on a reanalysis of the original data, per- 
formed in Ref. IsH . The original analysis contained an error in the 
kinematic treatment of the decay photon. See Ref. [3lj for details. 

''The values for Ml width are based on a reanalysis of the original 
data, performed in Ref. 3lj, using the values of <5i and Sq obtained 
experimentally in Ref. [31^1 . 



|31| which indicates a negligible second-forbidden contri- 
bution to the decay, in agreement with shell model pre- 
dictions. The early experimental determinations of Fmi 
^2f .^J disagree with the recent and most precise result 
|3ll | by about 2a. The recent result [31| is in best agree- 
ment with f3-a angular correlation experiments, as will 
be discussed later, and is adopted in this work. 



f{E,) = -SME.) 



9iE.) 



2 f2Amn 

En 



f{E.)- 



(12) 



(13) 



The isoscalar radiative widths are not related to /3 decay 
from factors by CVC. 



The matrix elements b(Ex) and c(Ex) have different 
functional dependences. This was first observed 28, 33 
through a comparison of the shapes of the final state dis- 
tributions in *Be following the a and 7 decays. The form 
of b(Ex) was later described JU using the R-matrix ap- 
proach, which parametrized ^(E^:) as an interfering sum 
of three different matrix elements, Aii, to the three 2+ 
levels in ^B shown in Fig. ^ 
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FIG. 8: The top panel shows the functional dependence of 
the weak magnetism matrix element, b(Ex), based on mea- 
surements of radiative decay in *Be. The bands indicate 
la experimental uncertainties. The bottom panel shows the 
Gamow- Teller matrix element, c(Ea;), based on fits to the al- 
pha spectrum discussed in Sec. Illll The uncertainties in c(Ex) 
are comparable to the width of the line and are negligible in 
the context of recoil order corrections. 



The notations used here are identical to those in Sec. 
mil We use the parameters reported in Ref. [Sj to de- 
termine 6(Ea;). The form of c(E^) was given in Eq. |S1 
and determined by fits to the a spectrum. We note that 
the R-matrix parameters appearing in both Eqs. Inland 
1141 may take different values in the two expressions. The 
forms of b(Ex) and c(Ex) are shown in Fig. |S1 



D. /3-q: Angular Correlations 

The (3-a angular correlations in the mirror decays of 
^Li and ^B have been measured several times as a func- 



tion of f3 particle energy [23, |32, |33j . Such measurements 
constrain the weak magnetism matrix element, 6, as well 
as the induced tensor, d. The angular correlations take 
the form 

N^[e,Ep,E,) ^l + a^{Ep,E,)cose+p^{Efi,E,)cos^9, 

(15) 
where the -(-I-) subscript refers to the *Li(*B) decay, 6 is 
the angle between the (3 and a particles, and the factor 
v/c for the f3 particle has been set equal to 1. The a^ 
coefficients are dominated by kinematic considerations, 
while the pip coefficients are strongly dependent on recoil 
order contributions. 
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where A=17.468 MeV is the total energy released in the 
^B f3-a decay chain. 

Assuming isospin symmetry, taking the sum and differ- 
ence of p~ and p+ produces cancellation between many 
of the mirror matrix elements of the ^B and ^Li decays. 



Corrections due to isospin breaking will be considered 
later. Defining d±= p_±p4_, dropping the vector matrix 
elements / and g, integrating over excitation energy E^; 
gives S± as a function of /3 particle energy, 



11 



S-{Ef3) 



Am,, ^ f b{E,)c{E^){A - E^ - EpfdE., 
Ep Jc^E,){A-E,-E0ydE, 



(17) 



s+iE^y- 



Am„ ^[c{E.)-d{E,) + ^^n{E.) 



2AmZ 



^^{E^) Ar, 



c{E,){A-E,-Ep)^dE, 



En 



Jc^iE,){A-E,-EpydE, 
I 



(18) 



12 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




' 


10 




[ ti ^ * --Li. 


" 1 1^ II ^U^ " 


CO. 

B 

< 4 
io 


1 


1 1 


1 h-HH 

ii Id I 


- 


2 


- 




- 





T 


, , , 1 . 


, , , . 



4 6 8 10 12 14 

Beta Particle Energy, Ep (MeV) 

FIG. 9: The solid squares are experimental data on 5~ from 
f3-a angular correlation measurements from Ref. 1331 . the 
open squares are from Ref. |23l . The curves indicate the 
la error bands from the prediction for 5^ based on Eq. 1171 
using the weak magnetism, bCEx), and Gamow- Teller, c{Ex), 
matrix elements. 



where the second class contribution to the induced tensor 
has been omitted. This is consistent with existing data in 
the A=8 nuclear system J31| , and with theoretical models 
which predict a second class current to contribute at a 
level below the current experimental sensitivity. 

The matrix elements &(Ej.) and c(Ej,), determined pre- 
viously, were applied to Eq. ll7l to predict the d~ observed 
in /3-a angular correlation measurements [23, Is^, 123 ■ 
The predictions are compared to the experimental S~ 
data graphically in Fig. El 

The level of agreement between the Eq. ll7lprediction. 
based on the radiative decay and alpha spectrum data, 
and the (3-a data sets was quantified by allowing the 
magnitude of b to float by a multiplicative constant, b — > 
«&, in Eq. 1171 The experimental data from the (3-a 
correlation measurements was then used to determine the 
best fit value of n. A value of k different than unity would 
indicate a disagreement between the radiative width data 
and the (3-a angular correlation data. This approach was 
previously applied 22, 28, 30, 31] with the motivation 
of testing CVC and searching for second-class currents. 
Here, the validity of CVC and the absence of second-class 
currents are assumed, and the test is performed to gauge 
the level of agreement between the two types of recoil 
order measurements. 



The best fit to the S angular correlation data from 
Ref. ^ gave k=1.06(4) with x-Vdof=7.7/8. The best 
fit to the data from Ref. 22] gave k=0.99(3) with 
X^/dof=24.8/15, where the large x^ value may be the 
result of the large point-to-point scatter of the data. 
The uncertainties in the data from Ref. [22 were ex- 
panded by \/x^/dof to account for this effect, and both 
data sets were fit simultaneously, yielding k=1. 014(26) 
with x^/dof=24.6/24. The values of k obtained, consis- 
tent with unity, indicate agreement between the radiative 
width measurement [3 111 and the (3-a angular correlation 
measurements [23, |33 , and provide confidence in the ex- 
tracted weak magnetism matrix element. 

The experimental i5^(E^) data is sensitive to the in- 
duced tensor matrix element, d. The effect of d on the 
neutrino spectrum is much milder than that of b. The 
energy dependences of b{Ejx) and c{Ex) were inferred 
directly from 7 and a spectrum measurements, respec- 
tively, but for d(Eix) there is no such experimental signal. 
The determination of the induced tensor is further com- 
plicated by the presence of the axial second-forbidden 
terms, J2 and J3, which appear in the expression for S~^ , 
Eq. 118! Fortunately, the infiuence of d on the neutrino 
spectrum is sufficiently small that very conservative esti- 
mates of uncertainty may be imposed on d without sig- 
nificantly inflating the total uncertainty of recoil order 
corrections. 

The (3 particle asymmetry from a polarized source of 
^Li or *B is also sensitive to J2 and J3, and would comple- 
ment P-a correlation measurements to allow a more pre- 
cise determination of the second- forbidden terms. One 
measurement of the asymmetry has been performed in 
®Li [33, but was systematically skewed by (3 particle 
scattering and required a sizable phenomenological cor- 
rection. We do not include the asymmetry measurement 
in our analysis, but note that future measurements of 
this type would be helpful in constraining the values of 
J2 and J3. 

Several models |23,|2J,|33 have been employed to esti- 
mate the magnitude of the axial second-forbidden terms. 
The models predict contributions to S^ from J2 and J3 
which are comparable to the contributions from the in- 
duced tensor, dj. It has been pointed out JS^] that 
mesonic exchange effects may be significant in A=8 /?- 
decays, especially at the second-forbidden level, and that 
shell model calculations may break down. 
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FIG. 10; The solid squares are experimental data on S^ 
from /3-Q angular correlation measurements from Ref. [33, 
the open squares are from Ref. |2^ . The curves indicate the 
la error bands from the prediction for 5+ based on Eq. 1171 
using the weak magnetism and Gamow- Teller matrix elements 
determined previously. Second forbidden contributions from 
j2 and J3 were ignored in this fit. 



To determine of the best value of d from the 5+ data, 
second-forbidden contributions are neglected and d will 
be assumed to take the same functional form as the 
Gamow- Teller matrix element, c. The possibility of large 
second- forbidden contributions to 6^, with magnitude 
given by the shell model predictions, will then be con- 
sidered and their effect on the extracted value of d will 
be assigned as an uncertainty. The uncertainty associ- 
ated with the ambiguity in the functional form of d will 
be estimated by fitting the S^ data with the assump- 
tion that d takes the same form as the weak magnetism 
operator, b. 

Utilizing the above assumptions, j2 and J3 are set equal 
to zero, d is considered to have the same form as c, d = 
rjc, and Eq. ^1 is used to fit the S~^ data, with 77 as 
the only parameter. The best fit to the (5+ data from 
Ref. ^31 gives ?7=10.3(2.3) with xVdof=2.7/8. The best 
fit to the data from Ref. 22] gives 77=10.6(1.4) with 
X^/dof=12.2/15. Fitting both data sets simultaneously 
gives 77=10.5(1.2) with xVdof=15.0/24. The results of 
the fits are compared to the J+ data in Fig. 1101 

The uncertainty associated with the second forbidden 
terms is estimated by assuming the values obtained using 



the model of Ref. 24], J2M 



400 and ijMc 



E. Recoil Order Effects on the Neutrino Spectrum 

The values and uncertainties of the weak magnetism, 
6, and induced tensor, d, matrix elements have been de- 
duced from experimental data. A further uncertainty 
is applied to these values due to imperfect isospin sym- 
metry and electromagnetic effects. The effect of isospin 
breaking is estimated by comparing the Gamow- Teller 
matrix elements of the ^B and *Li mirror /3 decays. Pre- 
vious comparisons of experimental a spectrum following 
*B and ^Li decays indicate Ch/cb ~1.07 ilS., ^]. As 
seen from Eqs. 1171 and 1181 this uncertainty propagates 
linearly to the extracted values of b and d. We thus as- 
sign to b(Ex) and d{Fix) a further 7% uncertainty, added 
in quadrature with previously stated uncertainties. Fur- 
ther electromagnetic effects, such as the difference in de- 
cay energies of ^Li and ^B and final state electromagnetic 
interactions, are discussed in Ref. |23 ] and are propor- 
tional to the second forbidden axial terms, J2 and J3. 
These effects contribute up to 4%, when the largest shell 
model values for j2 and j^ are assumed. We thus add, in 
quadrature, a further 4% uncertainty to b(Eix) and d{Ex)- 

Fig. Illlshows the ratio of b(Ex) to c(Ex) over the range 
of allowed excitation energies in ^Be. At high excitation 
energies, c{Ex) increases rapidly while b{Ex) decreases, 
as can be seen in Fig. |H1 In terms of the R-matrix 
approach, this is explained by comparing the Gamow- 
Teller strength of the high-lying doublet to the strength 
of the first excited state at 3.0 MeV, A^a/A^i=-11.8(8). 
For the weak magnetism transition, the ratio is much 
smaller, Ma/Mi=1A{1.6) [SJ, and the doublet transi- 
tion strength plays a smaller role. At excitation energies 
above 3.0 MeV, the result is a constructive interference 
of the Ma and the A^i terms for c(Ex)- Conversely, for 
excitation energies below 3 MeV, the terms interfere de- 
structively, causing c(Ex) to drop off more rapidly than 
b{Ex) and increasing the ratio b(Ej:)/c(Ex). 

The induced pseudoscalar matrix element may be esti- 
mated by applying the partially conserved axial current 
hypothesis, which indicates 



h{Ex) 



AM' 



mi, 



c{Ex). 



(19) 



- 750. We take d = rjc and the (5+ data |2S 
yielding 7/=13.8(1.2) with xVdof=16.5/24. 



are fit. 



The induced pseudoscalar appears only in the last term of 
1101 which is suppressed by a factor m^/M^. The induced 
pseudoscalar contribution to the /3 and neutrino energy 
spectra is thus of order m^/m^, and is ignored. 

The magnitude of recoil order effects on the *B neu- 
trino spectrum determined by the present treatment is 
compared to the previous treatment PJ] in Fig. 1121 



The uncertainty associated with the unknown func- 
tional form of the induced tensor is estimated by tak- 
ing d = ^b. A simultaneous fit to the (5+ data sets 
|22l I33, assuming no second forbidden contributions, 
gives e=0. 185(20) with xVdof=15.3/24. 



V. DETERMINATION OF THE *B NEUTRINO 
AND POSITRON SPECTRA 

Radiative corrections to nuclear /3 decay were first ex- 
plicitly formulated in [40j, and are exact to 0{a), where 
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FIG. 11: The solid curves indicate the la error bands in the 
ratio between the weak magnetism matrix element, fe(Ea;), 
and Gamow- Teller matrix element, cCEx), used in this work. 
The dashed line represents the ratio from Ref. JS^L used in 
previous determinations of the neutrino spectrum P.l9l l34ll3^ 
which neglected the excitation energy dependence of b(Ex) 
and c(Ea;). 
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FIG. 12: (Top panel) The normalized neutrino spectrum 
deduced in this work. (Bottom panel) The magnitude of the 
effect of the recoil order matrix elements on the neutrino spec- 
trum. The grey region shows the ztlcr band of the results 
obtained in this work. The black line was obtained using the 
recommended values from Ref. 34], which have been used in 
previous determinations of the neutrino spectrum. 



FIG. 13: (Top panel) Comparison of the direct positron spec- 
trum measurement 9\ with the predicted spectrum based on 
the measured alpha spectrum. The amplitude of the pre- 
dicted spectrum was floated. (Bottom panel) Same compari- 
son, showing only the momentum raiige measured in the di- 
rect positron spectrum measurement Q. 



The /?+ decay strength function determined in Sec. IIIII 
was applied, using Eq. 1 101 to determine the positron and 
neutrino spectra of ^B. 

The deduced positron spectrum was compared to 
the experimental spectrum |2j, and a one parame- 
ter fit to determine the amplitude gave an agreement 
of x^/dof=33.1/31, where only statistical uncertainties 
were included in the minimization function. The agree- 
ment is shown in Fig. 1131 The deduced positron spec- 
trum was then allowed to float by an energy offset and 
marginal improvement (x^/dof=32.6/31) was found for 
an offset of -14±20 keV. We note that the values quoted 
here are slightly different than those in Ref. Q due to the 
improved treatment of recoil order effects. The calibra- 
tion uncertainty of the positron measurement is reported 
as 25 keV |33, and is not included in the fits. 

The neutrino spectrum is presented numerically, with 
uncertainties, in Table Hvl 



a is the electromagnetic fine structure constant. Further 
corrections, dependent on the structure of the nucleus, 
occur at the O(a^ln^) level. These model dependent 
corrections are insignificant when compared to the exper- 
imental uncertainties in the neutrino spectrum and are 
not included. Radiative corrections for the case where 
the neutrino is detected while the positron remains un- 
observed were calculated explicitly in Ref. -41] , and affect 
the *B neutrino spectrum at the level of 1%. 



VI. CONCLUSIONS 

An accurate determination of the ®B neutrino spec- 
trum is important for the proper analysis of solar neu- 
trino data. Measurements of the a energy spectrum fol- 
lowing the *B /3+ decay provide the most direct method 
of inferring the /S"*" decay strength function, which is used 
to predict the neutrino spectrum. The a spectrum exper- 
iment described in this work was designed to eliminate 
several of the systematic effects common to past mea- 



TABLE IV: The neutrino spectrum of '^B 
in a given energy range. The spectrum is 
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and its uncertainties. Here P=dN/dE^ is the probability of a neutrino being emitted 
normalized to 1000 when integrated in terms of MeV. 



E. 



AP 



E. 



AP 



E. 



AP 



E. 



AP 



E, 



AP 



0.10 


0.214 


0.020 


3.20 


77.526 


0.299 


6.30 


130.963 


0.091 


9.40 


92.857 


-0.244 


12.50 


19.735 


-0.165 


0.20 


0.763 


0.043 


3.30 


80.456 


0.301 


6.40 


131.101 


0.074 


9.50 


90.528 


-0.248 


12.60 


17.963 


-0.157 


0.30 


1.513 


0.013 


3.40 


83.337 


0.301 


6.50 


131.134 


0.057 


9.60 


88.161 


-0.252 


12.70 


16.266 


-0.149 


0.40 


2.507 


0.021 


3.50 


86.164 


0.302 


6.60 


131.063 


0.040 


9.70 


85.759 


-0.255 


12.80 


14.647 


-0.140 


0.50 


3.763 


0.031 


3.60 


88.931 


0.301 


6.70 


130.888 


0.023 


9.80 


83.328 


-0.258 


12.90 


13.110 


-0.132 


0.60 


5.239 


0.041 


3.70 


91.635 


0.300 


6.80 


130.611 


0.007 


9.90 


80.869 


-0.261 


13.00 


11.655 


-0.123 


0.70 


6.914 


0.053 


3.80 


94.272 


0.298 


6.90 


130.232 


-0.010 


10.00 


78.387 


-0.263 


13.10 


10.286 


-0.115 


0.80 


8.772 


0.065 


3.90 


96.839 


0.296 


7.00 


129.752 


-0.027 


10.10 


75.885 


-0.264 


13.20 


9.005 


-0.106 


0.90 


10.798 


0.077 


4.00 


99.331 


0.292 


7.10 


129.174 


-0.039 


10.20 


73.368 


-0.265 


13.30 


7.813 


-0.097 


1.00 


12.976 


0.091 


4.10 


101.746 


0.288 


7.20 


128.497 


-0.051 


10.30 


70.837 


-0.265 


13.40 


6.712 


-0.088 


1.10 


15.292 


0.104 


4.20 


104.081 


0.284 


7.30 


127.724 


-0.063 


10.40 


68.298 


-0.265 


13.50 


5.703 


-0.080 


1.20 


17.735 


0.118 


4.30 


106.332 


0.278 


7.40 


126.856 


-0.075 


10.50 


65.754 


-0.264 


13.60 


4.787 


-0.071 


1.30 


20.292 


0.132 


4.40 


108.497 


0.273 


7.50 


125.895 


-0.087 


10.60 


63.209 


-0.263 


13.70 


3.965 


-0.062 


1.40 


22.950 


0.145 


4.50 


110.574 


0.266 


7.60 


124.843 


-0.097 


10.70 


60.667 


-0.262 


13.80 


3.237 


-0.054 


1.50 


25.699 


0.158 


4.60 


112.560 


0.259 


7.70 


123.701 


-0.108 


10.80 


58.131 


-0.259 


13.90 


2.602 


-0.046 


1.60 


28.528 


0.172 


4.70 


114.452 


0.252 


7.80 


122.471 


-0.118 


10.90 


55.606 


-0.257 


14.00 


2.058 


-0.038 


1.70 


31.427 


0.184 


4.80 


116.250 


0.244 


7.90 


121.156 


-0.128 


11.00 


53.095 


-0.254 


14.10 


1.602 


-0.031 


1.80 


34.386 


0.197 


4.90 


117.951 


0.236 


8.00 


119.758 


-0.138 


11.10 


50.602 


-0.251 


14.20 


1.228 


-0.024 


1.90 


37.395 


0.208 


5.00 


119.553 


0.227 


8.10 


118.278 


-0.148 


11.20 


48.131 


-0.247 


14.30 


0.929 


-0.019 


2.00 


40.447 


0.219 


5.10 


121.056 


0.218 


8.20 


116.720 


-0.158 


11.30 


45.686 


-0.242 


14.40 


0.694 


-0.014 


2.10 


43.531 


0.230 


5.20 


122.457 


0.209 


8.30 


115.086 


-0.166 


11.40 


43.271 


-0.238 


14.50 


0.513 


-0.011 


2.20 


46.640 


0.240 


5.30 


123.755 


0.199 


8.40 


113.378 


-0.175 


11.50 


40.889 


-0.233 


14.60 


0.376 


-0.008 


2.30 


49.767 


0.249 


5.40 


124.951 


0.189 


8.50 


111.599 


-0.185 


11.60 


38.545 


-0.227 


14.70 


0.273 


-0.006 


2.40 


52.903 


0.258 


5.50 


126.042 


0.179 


8.60 


109.751 


-0.193 


11.70 


36.242 


-0.222 


14.80 


0.196 


-0.004 


2.50 


56.041 


0.266 


5.60 


127.028 


0.168 


8.70 


107.838 


-0.200 


11.80 


33.984 


-0.215 


14.90 


0.140 


-0.003 


2.60 


59.174 


0.273 


5.70 


127.909 


0.157 


8.80 


105.862 


-0.208 


11.90 


31.774 


-0.209 


15.00 


0.099 


-0.002 


2.70 


62.296 


0.279 


5.80 


128.683 


0.146 


8.90 


103.827 


-0.215 


12.00 


29.616 


-0.202 


15.10 


0.069 


-0.002 


2.80 


65.401 


0.285 


5.90 


129.351 


0.135 


9.00 


101.734 


-0.222 


12.10 


27.515 


-0.195 


15.20 


0.047 


-0.001 


2.90 


68.482 


0.289 


6.00 


129.914 


0.124 


9.10 


99.587 


-0.228 


12.20 


25.472 


-0.188 


15.30 


0.000 


0.000 


3.00 


71.533 


0.293 


6.10 


130.369 


0.113 


9.20 


97.390 


-0.234 


12.30 


23.493 


-0.181 


15.40 


0.000 


0.000 


3.10 


74.549 


0.296 


6.20 


130.719 


0.102 


9.30 


95.146 


-0.239 


12.40 


21.579 


-0.173 


15.50 


0.000 


0.000 



surements, and precisely determines the strength func- 
tion characterized by the many level R-matrix approxi- 
mation. 

The primary uncertainty in the a spectrum measure- 
ment arises from the determination of energy scale, which 
was dominated by the temporal gain shift over the seven 
day run. Uncertainties in calibration from the implanted 
^°Na source and the correction for positron energy loss in 
the detector also contribute. Uncertainty in the energy 
scale near the spectrum peak was 9 keV. The uncertain- 
ties in the a spectrum were propagated to the neutrino 
spectrum using the R-matrix approach. 

Recoil order effects provide a significant contribution 
to the *B /3 decay, and have been treated in this work 
using the best available experimental data. The results 
differ from the previous treatment |3J| which has been ap- 
plied in recent determinations of the neutrino spectrum 



The primary component of the recoil order corrections 
is the weak magnetism term, b. Uncertainties in b are due 
to experimental effects l31| as well as isospin breaking 
and electromagnetic effects, and are included as uncer- 
tainties in the neutrino spectrum. Uncertainties in the 
neutrino spectrum from recoil order effects are roughly 
half as large as the uncertainties from the a spectrum 
measurement. 

The a spectrum reported here is in substantial dis- 
agreement with the previous measurement of comparable 
precision |g], but is in good agreement with the direct 
measurement of the positron spectrum |9j . 

This work was supported by the Department of Energy 
under Contract Nos. W-31-109-ENG-38 and DE-AC03- 
76SF00098. 
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